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Abstract 

In quantum state tomography, one potential source of error is uncontrolled contact of the system 
with a heat bath whose detailed properties are not known, and whose impact on the system more- 
over varies between different runs of the experiment. Precisely these variations provide a handle 
for reconstructing the system's effective relaxation dynamics. I propose a pertinent estimation 
scheme which is based on a steepest-descent ansatz and maximum likelihood. After reconstruct- 
ing the relaxation dynamics, the original quantum state of the system can be constrained to a 
curve in state space. 
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1. Motivation 

Quantum state tomography |1] is prone to errors, of various origins. For instance, samples 
might be small, the observables measured might not be informationally complete, or measure- 
ment devices might be inaccurate. In this letter, I focus on one further source of systematic 
error: The quantum system might be coupled in an uncontrollable fashion to a heat bath of fixed 
but unknown temperature (say, a piece of equipment or the surrounding air), whose effect may 
moreover vary randomly from sample to sample (say, because the duration of interaction varies 
between different runs of the experiment). Transient contact with the bath then triggers partial 
relaxation (to varying extent) of the system towards thermal equilibrium. In general, neither the 
original quantum state of the system nor the pertinent system-bath dynamics are known; and 
so neither is the target equilibrium state, let alone the trajectory in state space leading from the 
original state towards the latter Under such adverse circumstances, how can one hope to learn 
anything about the original quantum state? In this letter, I provide a simple, approximate scheme 
that allows one at least to constrain the unknown quantum state to a one-dimensional submanifold 
of state space, based on the outcomes of multiple independent runs of the experiment. The re- 
construction of this one-dimensional constraint manifold assumes a generic relaxation dynamics 
and uses the maximum likelihood approximation. 
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2. Generic relaxation dynamics 



In this section I first discuss the generic relaxation dynamics close to equilibrium and subse- 
quently propose an extension to states further away. Close to equilibrium, in the linear response 
regime, states typically relax towards equilibrium exponentially, 

p(f)-cr = exp(-f/Trel)[p(0)-Cr], (1) 

on some pertinent time scale Tiei. Here cr denotes the equilibrium state of the system. In the 
setting considered here, both Tiei and cr are fixed but unknown. The relaxation time Tiei should 
be larger than the typical contact times with the bath so that relaxation of the system towards 
equilibrium is only partial. 

Such exponential relaxation can be regarded as a steepest-descent flow of the relative entropy 
S{p\\cr), in the following sense. Given any set of observables {Ft,} which is informationally 
complete, every state p can be written in the Gibbs form 



p = Z(/l) ' exp 



(lncr-(ln(r)J-^/Fi 

b 



(2) 



with properly adjusted Lagrange parameters and partition function Z{A) 13-0] ■ (This is not 
to be confused with a thermal state: A thermal state has the Gibbs form, too, yet with only one 
observable (Hamiltonian) and Lagrange parameter (inverse temperature) in the exponent, and 
with cr replaced by the totally mixed state.) The Lagrange parameters [A^] and the expectation 
values [ft, :- {Fb)p}, respectively, constitute two possible choices for the coordinates in state 
space. Associated with these coordinates are the respective basis vectors [dh ■- d/dA''} and 
:= d/dfb}, as well as their one-form duals [dA"} and {dfa}, related via dA"{dh) - dfbid") - 6"^ 
list]. The transformation from one coordinate basis to the other is effected by the correlation 
matrix Cab d^\nZ{A)ldA"dA\ 

da = -Cabd' , d" = -iC-'ydt (3) 



(where I have adopted the convention that identical upper and lower indices are to be summed 
over), and likewise the transformation of their duals, 

dfa = -CabdA' , dA" = -(C-'r'dfb. (4) 

Being symmetric and positive, the correlation matrix also provides a natural Riemannian metric 
on state space, the Bogoliubov-Kubo-Mori (BKM) metric |6-11]. Its elements may be regarded 
as the components of the metric tensor 

C := C„it/r ®(//. (5) 

The relative entropy S(p\\cr) has the gradient 101 

dS(p\\cr)^-A"dfa, (6) 

and so the pertinent steepest-descent curve (in the BKM metric introduced above) has a tangent 
vector 

V^C-H-dS(p\\cr))^-A"da. (7) 
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Consequently, on a steepest-descent trajectory the Lagrange parameters must evolve according to 
dA'' jdt oc {-A"da)A'' - -/t* and thus relax exponentially towards A-Q, which corresponds to the 
state cr. In the linear response regime, exponential relaxation of the Lagrange parameters in turn 
implies exponential relaxation of expectation values, and so indeed of the statistical operator, 
Q.E.D. 

Identifying thus the generic relaxation dynamics near equilibrium with a steepest-descent 
flow opens the possibility to formulate a generic dynamics even in regions further away: In the 
absence of detailed information about the system-bath dynamics, and regardless of how far the 
state of the system might be from equilibrium, I henceforth assume that the generic effect of a 
bath on the system is to move its state by some (unknown) distance along the steepest-descent tra- 
jectory towards some (unknown) equilibrium state cr. Such an extension of the steepest-descent 
paradigm to full state space has been studied before, and in the special case of the classical 
Fokker-Planck equation (with a particular choice of metric) has been shown to reproduce the ef- 



fective dynamics excactly I121I13I1 . In the setting considered here, the steepest-descent algorithm 
induces an effective quantum operation on the state space of the system. Close to equilibrium, 
this effective quantum operation is a completely positive map. Further away from equilibrium, 
however, it might no longer be linear and thus not completely positive; which suggests that, in 
contrast to a completely positive map, the steepest-descent ansatz may presume non-negligible 
initial correlations between system and bath. Such initial correlations are not surprising in view 
of the fact that the properties of the heat bath are not known a priori; they simply reflect the 
fact that learning something about the initial state of the system will entail learning something 
about the state of the bath, too. (There are also other reasons why one should expect the effective 
quantum dynamics to be nonlinear, see e.g. Ref. 1 14].) It remains to be seen whether for a given 
steepest-descent flow there exists always some compatible microscopic system-bath dynamics, 
just as any completely positive map always arises from some underlying unitary evolution of the 
(initially uncorrected) composite system-bath state. I conjecture that this is the case, but defer 
its proof to future work. 

Any steepest-descent trajectory towards cr can be parametrized in the form 

p(y) oc exp [(In cr - (In cr>^) - jG] , (8) 

with some single observable G and associated Lagrange parameter y. I prove this assertion 
in two steps: (i) any curve parametrized in this form is a steepest-descent curve towards cr; 
and (ii) for every initial state po there exists a curve of this form passing through it. The first 
statement becomes evident when one considers the above parametric form to be a special case 
of the more general form (|2), with the pair (y, G) being one of the informationally complete 
{{A'', Fb)} and the Lagrange parameters associated with all Ft, G equal to zero. If these other 
Lagrange parameters are initially zero, they remain zero under the steepest-descent dynamics 
dA''/dt oc -A''. So indeed, a curve parametrized as above describes a steepest descent towards cr. 
As for the second statement, for any initial po (except for states on the boundary of state space) 
the choice G = In cr - In po will yield a steepest-descent curve passing through both po (at y = 1) 
and cr (y = 0). Given po and cr, this choice for G is unique up to multiplicative and additive 
constants. 

Given the imperfect data from multiple runs of the experiment that have been affected (to 
varying degree) by uncontrolled contact of the system with the same unknown heat bath, a precise 
reconstruction of the original quantum state is clearly not possible. However, assuming that 
transient contact with the bath triggers a generic relaxation dynamics as described above, it is at 
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least possible to constrain the original quantum state to some steepest-descent trajectory. Having 
established the parametric form (|8) of such a trajectory, the task is then to infer from the available 
experimental data the observable G (up to multiplicative and additive constants) as well as some 
state (not necessarily the target state cr) which lies on this trajectory. This is the inference task to 
which I turn in the next section. 



3. Inferring the constraint curve 

The experiment is run multiple times. In the /th run, after the supposed disturbance by the 
bath, one performs complete quantum state tomography on a sample of size A^, and obtains the 
tomographic image ju,. Different runs of the experiment are disturbed by the same bath but to 
varying extent (say, because the contact time varies between runs). As a result, the original 
quantum state is shifted on its relaxation trajectory towards equilibrium to varying degree; and 
so rather than being clustered around a single point (which would yield a unique state estimate), 
the images are expected to be spread out along this trajectory. 

Assuming that the relaxation trajectory has the steepest-descent form (|8]l, the problem of 
inferring the observable G from images that correspond to various values of y is analogous to the 
problem of estimating a Hamiltonian from thermal data at various temperatures. Assuming that 
the tomographic images are not too far apart, and expanding G in terms of the informationally 
complete set {Fb}, 

G = -^Ft (9) 

(with the same index summation convention as above), the log-likelihood of observing the data 
{/i,) is asymptotically (i.e., for sufficiently large sample sizes) given by 

L({//,)|f, a) ~ (A^/2)[(r)^ - {C-'Y'6fa{^, cr)6M^, cr)], (10) 



modulo additive terms that do not depend on ^ or cr u3 . Here := 2, and 



(D, (11) 
is the "expectation value" of the covariance matrix 

The latter in turn depends on the deviations of the sample means :- {Fi,)^. from their weighted 
averages ft 2; '^ifl^ with the rth sample carrying relative weight w, := A^,7A^. The correlation 
matrix C„;, is evaluated at the center of mass Jl :- 2, w,yu, of the tomographic images. Finally, 
the valuations 6f are given by 

6m,cT):^{Fh)n^^,^)-fb, (13) 

where ni^, cr) is the unique state of the form ([8]l which yields (G),r = (G>p. 

Unlike in the problem of estimating a Hamiltonian from thermal data where cr corresponds 
to a given reference state, the target state cr is not fixed a priori but is itself a variable to be 
inferred. To maximize the above log-likelihood as a function of both ^ and cr, the steepest-descent 
trajectory must maximize (F)^ as well as render 6f -0. The first requirement implies that ^ must 
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be the dominant eigenvector (i.e., the eigenvector associated with the largest eigenvalue) of C"'r, 
the product of the inverse correlation and covariance matrices. Such an eigenvector condition is 
a familiar result in principal component analysis The second requirement, on the other 

hand, is met by any cr that ensures n - fi, and hence by any steepest-descent trajectory that 
passes through the center of mass fi. Thus one can infer from the experimental data both a 
maximum likelihood estimate for ^, and hence for G, and a state (namely ft) through which the 
steepest-descent trajectory must pass. Together, these characterize the steepest-descent trajectory 
uniquely. The original quantum state of the system must lie somewhere on this trajectory. 



4. Example: qubits 

The simplest example is state tomography on an exchangeable sequence of qubits, emitted 
by some i.i.d. quantum source. The tomographic experiment is performed several times, with the 
ith run yielding the tomographic image /i,-. Rather than isotropically around a point, these images 
[fii] turn out to be clustered along a curve in state space; say, along a straight line parallel to the 
y axis of the Bloch sphere (Fig. [TJ. This suggests that in between runs of the experiment, there is 
some uncontrolled change in the environment which is characterized by a single parameter An 
obvious candidate is the coupling to an unknown heat bath, with the contact time fluctuating in 
between runs. 

For a qubit, the informationally complete set of observables {Fi,} can be taken to be the Pauli 
matrices. The associated parameter vector ^ points in the direction of the effective "dissipative 
force" which drives the qubit towards equilibrium. Due to the nontrivial geometry of quantum 
states this force is in general not parallel to the spatial orientation of the data cluster (here: the y 
axis) as one might naively expect. Rather, it is tilted against this axis by some angle 0. This angle 
can be calculated as follows. Assuming that the measured data are perfectly aligned parallel to 
the y axis as shown in Fig. [T] all entries of the covariance matrix vanish except for Fyy 0. And 
for the center-of-mass state ji located in the x-y plane as shown in Fig. [1] the inverse correlation 
matrix has the form 



c-} C„i 
C:} 



(14) 



The requirement that f be the dominant eigenvector of C then implies f ; = and ^xl^y - 
C^,f /C,Ty . As a consequence, f is tilted against the y axis by the angle 

= arctan(C;,'/C^'). (15) 

This tilting angle vanishes whenever the center of mass lies on one of the two (x or y) axes. 
Away from the axes, however, the tilting angle is non-zero, and increases as the center of mass 
moves closer to the surface of the Bloch sphere. To illustrate the latter, I consider a center-of- 
mass state on the x-y diagonal, with azimuth 7r/4 and variable Bloch vector length r. The 
relevant elements of the inverse correlation matrix are then C,7y - 1/(1 - r^) and C^^f = Cy^ — 
(tanh"' r)/r. This yields a function 0(r) which increases monotonically from 0(0) = to 0(1) = 
7r/4, and which to a reasonable degree of accuracy can be approximated by 0(r) « {nlA)p-. 
Only near the origin of the Bloch sphere, therefore, and hence for highly mixed states, does the 
inferred direction of the effective dissipative force coincide with the "naive" estimate based on 
the spatial orientation of the data cluster. For states that are (nearly) pure, on the other hand, the 
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Figure 1: First quadrant of a two-dimensional section (z = 0) of the Bloch sphere. The black dots indicate the tomo- 
graphic images obtained in different runs of the experiment, and the smaU circle their center of mass fi. The Bloch 
vector associated with the center of mass is assimied to have no z component ((cr^)^ = 0), azimuth njA, and length r. 
Even though the images are aligned parallel to the y axis, the inferred direction f of the effective force which drives the 
relaxation towards equilibrium is tilted against that axis by an angle <p. This angle grows from (^ = Oatr = Oto0 = ;r/4 
at r = 1. 

reconstruction scheme presented here may yield a direction ^ which differs significantly from 
that naive estimate. 

5. Conclusions 

Whenever an i.i.d. sequence of quantum systems is subjected to disturbance by the same 
heat bath, yet to an extent that varies randomly between samples, multiple runs of quantum 
state tomography are expected to yield data spread out along a curve in state space. From the 
shape and location of this sprawl, one can infer the system's effective relaxation dynamics under 
the influence of the bath. Owing to the nontrivial geometry of quantum states, the result of 
this inference can be rather counterintuitive. In the qubit example, the inferred direction of the 
effective dissipative force generally deviated from the principal axis of the data cluster 

There are several ways to extend the results of the present letter, which are left to future 
work, (i) So far the reconstruction scheme for the effective relaxation dynamics does not take 
into account the prior distribution prob(^, cr) of its parameters. Doing so will lead to a Bayesian 
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modification of the maximum likelihood framework presented here. A full Bayesian analysis 
should also include a study of the error bars on ^ and cr. (ii) Many physical systems, especially 
larger ones, exhibit not a single relaxation time but a whole hierarchy of time scales pertaining to 
the relaxation of different sets of degrees of freedom. In this case thermalization occurs in stages, 
on successively longer time scales I2lll22ll . Provided the typical contact times with the bath are 
shorter than the smallest scale in that hierarchy, the reconstruction scheme presented here is still 
valid. It pertains then to the first thermalization stage, with cr being no longer the equilibrium 
state but the target state of this first stage. If contact times vary widely across different relaxation 
time scales, however, the scheme must be adapted, (iii) The system might be disturbed by not 
just one but several baths, each to independently varying degree, or the impact of a single bath 
might be governed by more than one parameter Then, rather than along a curve, tomographic 
images will be spread out in some higher-dimensional submanifold of state space. Inferring the 
dimensionality of this submanifold, as well as the associated parameters, will require a further 
generalisation of the scheme presented here, (iv) Reconstructing the relaxation trajectory (or sub- 
manifold) is the first step towards the reconstruction of the system's unperturbed quantum state. 
I argued that this unperturbed quantum state must lie somewhere on the relaxation trajectory (or 
submanifold); but its precise location will depend on further assumptions, not considered in the 
present letter, about the distribution of the contact times and possibly other parameters of the 
bath, (v) Finally, I consider it worthwhile to study in more detail the conceptual underpinning 
as well as possible modifications of the steepest-descent paradigm. And on a speculative note, 
taking the steepest-descent paradigm at face value and identifying (some suitable function of) 
relative entropy as actual "time," one might even be tempted to try to establish a link to (equally 
speculative) ideas in other areas of physics that aim to reduce the notion of time to a distance 
between configurations 1I23I1 or thermal properties [241. 
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